Granular Avalanches in Fluids 
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Three regimes of granular avalanches in fluids are put in light depending on the Stokes number 
St which prescribes the relative importance of grain inertia and fluid viscous effects, and on the 
grain/fluid density ratio r. In gas (r ^ 1 and St > 1, e.g., the dry case), the amplitude and time 
duration of avalanches do not depend on any fluid effect. In liquids (r ~ 1), for decreasing St, 
the amplitude decreases and the time duration increases, exploring an inertial regime and a viscous 
regime. These regimes are described by the analysis of the elementary motion of one grain. 

PACS numbers: 45.70.-n, 45.70.Ht, 46.10.+z,47.55.Kf. 
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Granular matter has received much attention from 
physicists over the past few years [Q . Beyond the funda- 
mental interest in the physics of granular systems which 
can present some features of either solids, liquids or even 
gases, the understanding of granular materials is essential 
in many industrial activities such as pharmacology, chem- 
ical engineering, food, agriculture, and so on. Many stud- 
ies concern the avalanches that may arise on the slope of 
a granular pile in air. Such granular avalanches occur 
in various places in Nature, from small scale, as for the 
building of any sand pile, to large scale, as the event ob- 
served after the Mont St-Helen eruption in 1980. Two 
angles can be defined when building a pile: the maxi- 
mum angle of stability at which an avalanche starts 
and the angle of repose 0r at which the avalanche stops. 
Between these two angles is a region of bistability where 
the grains can either be flowing ("liquid state") or at 
rest ("solid state"). Many experiments performed with 
dry grains in a rotating cylinder ^, ^, |[ ^ showed 
clearly the existence of these two angles. 

To date, no detailed study has focused on the influ- 
ence of the interstitial fluid for a totally immersed grain 
assembly. This influence is certainly important in granu- 
lar avalanche processes, as evidenced by the marked dif- 
ferences observed by geologists between subaqueous and 
eolian cross strata j7|. As a matter of fact, the propa- 
gation of subaqueous dunes differs in general from the 
propagation of eolian dunes even if the slope angles are 
quite similar: When the transport rate of sand particles 
is large enough, the flow is continuous in the lee side of 
the structure in the immersed case, but occurs by suc- 
cessive avalanches in the dry case Q]. This observation 
prompted geologists to accumulate data on avalanches of 
sand or beads in rotating drums filled with air or wa- 
ter Q or even with glycerol mixtures Q , that seemed to 
show that the amplitude of avalanches decreases and the 
time duration increases with the fluid viscosity. We have 
performed an extensive series of experiments to investi- 



gate the influence of the interstitial fluid on the packing 
stability and the avalanche dynamics. The analysis of 
our results obtained with a rotating drum set-up indi- 
cate the existence of three regimes: (i) a free-fall regime 
for which there is no fluid influence and that corresponds 
to the classical dry regime, and two regimes where the 
interstitial fluid governs the avalanche dynamics, namely 
(ii) a viscous regime and (iii) an inertial regime. 

Our rotating drum consists of a cylinder of inner di- 
ameter D ranging from 8 cm to 46 cm and lying on two 
parallel rotation axes. It is driven by a microstep motor 
followed by a 1/100 reducer and a rubber transmission 
so that the cylinder turns at the rotation rate Q by step 
of lO"'^ degree without shocks. The cylinder is half filled 
with sieved solid spheres of diameter d and density ps, 
totally immersed in a fiuid of density pf and dynamic 
viscosity fj,. The pile is confined between two parallel 
glass endwalls separated by the gap width b. In each 
experimental configuration, we used a sufficiently large 
gap width {b/d > 15). This study was achieved at low 
enough rotation rate fl to be in the intermittent regime 
of macroscopic avalanches |Q, ^ : The pile slope increases 
linearly with time at the rate fi, and then quickly relaxes 
by a surface avalanche process. A CCD camera aligned 
along the axis of the cylinder allows for visualization of 
the rotating pile. Images are taken at regular time inter- 
vals and then analyzed to track the pile interface. Except 
during the avalanches, the interface is found to be linear, 
with a roughness of the order of one grain size, which 
justifies the calculation of the averaged slope angle 6 of 
the pile with the resolution 0.01°. 

Two typical recordings for the time evolution of the 
average slope angle 9 are displayed in Fig. 1 for 230 
^m glass beads immersed either in air or water. We fo- 
cus on two typical parameters of the avalanche dynamics 
that we found uncorrelated: The avalanche amplitude 
characterized by the hysteresis angle A9 = 9,^ - 9r, and 
the avalanche time duration T which is calculated as the 
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FIG. 1: Time evolution of tlie average slope angle of a pile 
of glass beads (d = 230 ± 30 fim) immersed in air (a) or water 
(b) and contained in a rotating drum. The insert is a zoom 
of an avalanche in the water case. 



time interval between 5% and 95% of the corresponding 
avalanche amplitude. The mean values A9 and T are 
then calculated for one experiment over all the succes- 
sive macroscopic avalanches that affect the entire slope. 
In the following, we will drop the mean bar notations for 
simplicity. Despite the dispersion, two different behaviors 
clearly appear: In the air case (Fig. la), the avalanche 
amplitude A9 is large (few degrees) and the avalanche 
duration T is small (typically one second), whereas in 
the water case (Fig. lb) A9 is small (less than one de- 
gree) and T is large (typically one minute). One crucial 
parameter of the phenomenon appears to be the particle 
diameter d when one looks in Fig. 2 at the evolution of T 
and A9 as a function of d when the grains are immersed 
either in air or in water. In the air case, there is no signif- 
icant dependency of T and A9 on d. By contrast, in the 
water case, the avalanche duration which is close to the 
air case for the larger d, increases first sligthly (T cx d^^^) 
then drastically (T oc d^) from typically 2 to 100 s when 
the grain size is decreased from 1 mm to 0.18 mm. Paral- 
lely, the avalanche amplitude A9 in water corresponds to 
the air value for the larger d but decreases close to zero 
when the grain size is decreased to 0.18 mm. 

The pertinent dimensionless parameters governing the 
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FIG. 2: Time duration T and amplitude A9 of macroscopic 
avalanches as a function of the grain diameter d for glass 
beads immersed in air (■) or water (•) in a rotating cylinder 
of diameter D = 16 cm. The error bars correspond to the 
standard deviation. 



avalanche dynamics in fluids can be inferred by consider- 
ing the elementary falling process of one solid grain on its 
neighbour from below in a fluid and under the action of 
gravity. Let us write the simplified following equation of 
motion for the grain of velocity u down the slope between 
two collisions: 



TrPsd-^-rr = -Apgdhin9 - 
6 dt 



(1) 



Starting from zero velocity, the grain increases its grain 
momentum at the rate (Tr/6)psd^{du/dt) under the ac- 
tion of its apparent weigth {iT/6)Apgd'^sm9 (where Ap 
= ps - Pf) minus a fluid drag force Fd. At this ba- 
sic stage, the solid friction force can be modelled by a 
dynamical Coulombic term which just reduces the ap- 
parent gravity. Two behaviors can be discussed accord- 
ing to the value of the particle Reynolds number: At 
Re ^ 1, is equal to the viscous Stokes force Sndriu 
so the grain will possibly reach its viscous limit velocity 
Uoov — Apgd^siii9/18r] in the characteristic time Tcv = 
Psd^ /ISr], i.e. for a characteristic distance Scv = TcvUoov 
At Re ^ 1, Fd is the inertial fluid force Cd{TT/6)d'^ pfu'^ , 
and the inertial characteristic time and distance are Td 
= (ps/p/)^/^(2psd/Ap5sin6')i/2 and = TdUooi where 
Uooi = {2Apgds\n.9 / pfY^^ is the inertial limit velocity 
(for simplicity we take here the drag coefficient as con- 
stant : Cd — I/tt ~ 0.3). By comparing the two charac- 
teristic distances 5cv and Sd with the elementary distance 
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beween two successive collisions taken as the grain diam- 
eter d, we introduce two dimensionless numbers, St and 
r, which govern the grain dynamics in this elementary 
falling process: 5cv/d = 2(Tct,/ryj)^ — 2St^ and Sd/d 
= 2{Tci/Tfff = 2r2 where r// = (2psc;/Apgsin6')i/2 jg 
the typical timescale of free falling of a grain over d, St 
= (l/18A/2)/9y^(Ap5sin6')i/2d3/2/r7 is the Stokes number 
which prescribes the relative importance of grain inertia 
and fluid viscous effects of and r = {ps/ PfY^'^ is related 
to the density ratio. Note that the Reynolds number cor- 
responds to the third time ratio, so that Re = TcvItci = 
St/r. For St 1 and r S> 1, the grain does not reach any 
limit regime: This is the "free-fall regime" . For St <C 1 
and r ^ 1 the grain reaches its limiting Stokes velocity: 
This is the "viscous limit regime" . For St ^ 1 and r <C 
1, the grain reaches its limiting inertial velocity: this is 
the "inertial limit regime". For St <C 1 and r ^ 1, the 
grain reaches one limit velocity depending on the Re = 
St/r value: for Re ^ 1 (resp. Re » 1) the hmit regime 
is the viscous (resp. inertial) one. The exact boundaries 
between the three domains in the (St, r) plane of Fig. 3 
will be precised further. In this diagram are reported all 
our data and other data j|, ^ , corresponding to different 
sphere materials (glass and Nylon) in different fluids (air, 
water, silicone oils or glycerol mixtures of different vis- 
cosities). All experimental results correspond to roughly 
two data lines in this diagram: one data line for the liq- 
uid case where r ~ 1 and St ranges from 0.2 to 40 and 
another data line for the air case where r ^ 40 and St 
ranges from 30 to 10^. 

Let us now look if the complex dynamics of macro- 
scopic granular avalanches in fluids can be related to 
these elementary falling processes. In the two limit 
regimes, one may reasonably suppose that the time du- 
ration T of a macroscopic granular avalanche will scale 
as D I d elementary fallings each of time duration djUaav, 
so that Tu — D/U oov = ISrjD / Apg(Psm9 in the viscous 
limit regime and T, = D/Uooi = Dp}^^/(2Ap5dsin6')i/2 
in the inertial limit regime. By the way, the two scalings 
T cx c?~2 and T (x d~'^/'^ observed in Fig. 2 for small and 
large d respectively correspond to these two predicted 
scalings. We have plotted T/T^ and T/Ti in Fig. 4a 
and 4b respectively and we observe that all data collapse 
onto the plateau T ~ 4r„ for St < 5 (Fig 4a), and the 
plateau T ~ 2Ti for St > 3 (Fig 4b). A careful analysis of 
each data point of Fig 4a, b allow us to identify the crit- 
ical Reynolds number Rcc — 2.5 for the viscous/inertial 
transition, and thus to draw the corresponding boundary 
line of slope 1 between the two corresponding domains in 
the log-log plane (St, r) of Fig. 3. 

In the elementary free-fall regime corresponding to the 
air case (St ^ 1 and r ^ 1), the same approach would 
lead to the avalanche time duration T = {D/d)Tff = 
D{2ps/Apgdam9y^'^ , which is not consistent with the 
non d dependence observed in Fig. 2. In addition, we 
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FIG. 3: Diagram of the different regimes for the elementary 
falling process of one grain in the (St, r) plane. The three 
regimes are separated by the three boundaries Stc = 10, Rec 
= 2.5 and = 4. Our results (filled symbols) correspond 
to glass beads in air (■) or water (•), or in silicone oils of 
different viscosities (a), and to Nylon spheres (♦) in air, water 
and silicone oil. Other results (open symbols) of Evesque [5] 
(□) and Allen [10] (o) correspond to glass beads in air and 
different water-glycerol mixtures. 



have also observed that T scales as D^/^ in the air case 
rather than as D. All these observations lead us to con- 
sider a dry avalanche as a global accelerated rush over 
D of macroscopic timescale Tfj = {2psD/ ApgamOy/"^ 
rather than a succession of elementary falling processes. 
This scaling is indeed observed as all data in air collapse 
onto the plateau T ~ 3r// for St > 30 (Fig. 4c). In all 
three regimes, the plateau values of Fig. 4a,b,c are not far 
from one, meaning that our crude approach catches the 
essential of the avalanche dynamics in fluids. Note that 
for much larger D inaccessible to laboratory experiments, 
such a free-fall regime will possibly not be observed as the 
inertial fluid force or a solid friction force will come into 
play at large grain velocity. 

Let us focus now on the avalanche amplitude A6. 
When plotting A9 as a function of the Stokes number 
St (Fig. 5), we obtain a single master curve. At large St 
(St > 20) AO has the constant plateau value AO ~ 3°±1° 
whereas A9 clearly decreases with St at small St (St < 
20). We also observe that the decrease of A9 at small 
St is essentially due to a decrease of the maximum an- 
gle of stability 0^, the angle of repose Or being almost 
unchanged (cf. Fig. 1). This can be explained by con- 
sidering the collision process beween immersed grains. 
Indeed, the Stokes number is the only parameter that 
governs the coefficient of restitution e for the immersed 
binary collision between solid grains JTot . Under the crit- 
ical value Stc — 10, the coefficient of restitution is zero: 
the collision is totally inelastic as all the kinetic energy 
of the grain is dissipated by the fluid during the collision 
process Above Stc, the coefficient of restitution in- 
creases quickly with St and becomes close to its maximal 
"dry" value Cdry ^ 1 for St > 100 Q. It is worth- 
noting that the curve e = f(St) is the same whatever 
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FIG. 4: Time duration T of macroscopic avalanches normal- 
ized either by (a) the viscous timescale Tv (b) the inertial 
timescale Ti or (c) the free-fall timescale T// as a function 
of the Stokes number St for different grains in different fluids 
and rotating cylinders of different diameter D. Same symbols 
as in Fig. 3. 



the density ratio r pO[ . Considering again the immersed 
granular avalanches in fluids, the grain kinetic energy will 
be totally dissipated by the fluid in the collision process 
at low St, with an all the more smooth collision when 
St evolves towards zero. The obtained packing is thus 
certainly all the more loose, i.e. with a lower packing 
fraction As the maximum angle of stability of a 

granular pile depends largely on the arrangement of the 
packing, decreasing with the packing fraction p^ , this 
explains the decrease of with St at low St. 

In addition, if one consider the critical Stokes value 
Stc ~ 10 independent of r as the boundary line between 
the accelerated regime and viscous limit regime in the 
(St, r) plane of Fig. 3, this leads to the critical density 
ratio Tc ~ 4 separating the free-fall regime and the in- 
ertial limit regime. As r values larger than A{ps/ pf > 
16) can hardly be reached experimentally for solid/liquid 
system, the free-fall regime corresponds only to solid/gas 
systems like, e.g., the dry granular avalanches. 

Finally, for large St (St > 20), all the events corre- 
spond to macroscopic avalanches that affect the entire 
slope and no small event are observed between two suc- 
cessive macroscopic avalanches: The size distribution is 
a gaussian curve centered on the value ~ 3°. This 
kind of distribution, classically found for dry granular 



<l 




FIG. 5: Amplitude AS of macroscopic avalanches of glass 
beads in different fluids as a function of the Stokes number 
St. Same symbols as in Fig 3. 



avalanches ^, ^, |[ |j, is incommensurate with the 
ideas of Self Organized Criticality (SOC) developed by 
Bak et al. [ p^ which would predict a power law distri- 
bution without any typical scale. The reason may be 
the dissipation rate of the system [16]: In the cellular 
automata models, which illustrate nicely the SOC, au- 
tomata are strongly overdamped whereas dry granular 
experiments are weakly dissipative. When introducing 
inertia in cellular automata models, a complex distribu- 
tion, mixing power law distribution for the small events 
and gaussian distribution for the large events, is obtained 
|]l5| . For decreasing St (St < 20), we observe experi- 
mentally together with the gaussian distribution of large 
events the appearance of numerous small events (affect- 
ing not all the slope). We expect such a behavior in the 
viscous regime of low St as the hysteresis A6' of the sys- 
tem goes to zero with St, which is a condition for the 
system to evolve towards criticality. But up to now, we 
have not enough resolution to characterize the size distri- 
bution of these small events, and to conclude if it obeys 
or not the power law related to SOC. In addition, the 
regime of intermittent avalanches is hard to obtain when 
the regime is more and more viscous as the time duration 
of avalanches diverges. 

By conclusion, we have shown the existence of three 
regimes (free-fall, inertial limit, and viscous limit) for 
granular avalanches in fluids, controlled by the Stokes 
number which measures the ratio of particle inertia to 
viscous fluid effects, and the density ratio. The time du- 
ration of the macroscopic avalanches that affect the en- 
tire slope have been predicted in all these three regimes. 
The amplitude of these macroscopic avalanches has been 
shown to be constant at high St while decreasing with St 
at low St. Finally, more refined experiments remain to be 
done to see if the system evolves towards criticality when 
St tends towards 0, i.e. for highly dissipative systems. 

We acknowledge B. Andreotti, S. Douady, D. Lhuillier, 
O. Pouliquen, E.J. Hinch, and G.M. Homsy for fruitful 
discussions. 
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